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Band calculations for Ce compounds with the AuCu3-type crystal structure were carried out 
on the basis of dynamical mean field theory (DMFT). The auxiliary impurity problem was 
solved by a method named NCA/ 2 vc (noncrossing approximation including the f 2 state as a 
vertex correction). The calculations take into account the crystal- field splitting, the spin-orbit 
interaction, and the correct exchange process of the f 1 — > f°,f 2 virtual excitation. These are 
necessary features in the quantitative band theory for Ce compounds and in the calculation 
of their excitation spectra. The results of applying the calculation to CePd3 and CeRli3 are 
presented as the first in a series of papers. The experimental results of the photoemission spec- 
trum (PES), the inverse PES, the angle-resolved PES, and the magnetic excitation spectra were 
reasonably reproduced by the first-principles DMFT band calculation. At low temperatures, 
the Fermi surface (FS) structure of CePd3 is similar to that of the band obtained by the local 
density approximation. It gradually changes into a form that is similar to the FS of LaPd3 as 
the temperature increases, since the 4/ band shifts to the high-energy side and the lifetime 
broadening becomes large. 

KEYWORDS: dynamical mean field theory, band theory, ARPES, magnetic excitation, CePd3, CeRfi3, 
Cauchy integral 
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1. Introduction 

Nonempirical band calculations for strongly correlated 
electron systems have been extensively developed on the 
basis of dynamical mean field theory (DMFT). 1 ' 2 ) The 
4/ electrons in Ce compounds are typical strongly corre- 
lated electrons. 3 ~ 5 ' Recently, a DMFT band calculation 
scheme for Ce compounds was developed in refs. 6 and 7. 
In the present work, it is applied to Ce compounds with 
the AuCu3-type crystal structure, which show a wide va- 
riety of 4/ electronic states from the most itinerant limit 
to the localized limit. 

The 4/ state splits into the j = 5/2 ground multi- 
plet and the j = 7/2 excited multiplet with a separa- 
tion of about 0.3 eV owing to the spin-orbit interaction 
(SOI). The multiplet shows crystal-field splitting (CFS) 
of the order of 100 K. In cubic crystals, the j = 5/2 
multiplet splits into the (j = 5/2)T7 doublet and the 
(j = 5/2)r§ quartet. Hereafter, we call them the Tj and 
T$ states, respectively. It is important to take account 
of the SOI and CFS effects in 4/ systems. 3) In DMFT, 
the correlated band electron problem is mapped onto the 
calculation of the single-particle excitation spectrum of 
the auxiliary impurity Anderson model in an effective 
medium. Reliable methods of solving the impurity An- 
derson model with CFS and the SOI effect are needed 
in the DMFT band calculation for 4/ compounds. A 
theory named NCA/ 2 vc (noncrossing approximation in- 
cluding the / 2 state as a vertex correction) has been 
developed 8-11 -* and combined with the linear muffin-tin 
orbital (LMTO) method 12 ' 13 ' to carry out the DMFT 
band calculation. 7 ' NCA/ 2 vc can include CFS and the 
SOI effect, and also the correct exchange process of the 
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/ — > f°,f 2 virtual excitation. The calculation gives an 
accurate order of the Kondo temperature (Tk). 

The DMFT band calculation will be applied in a series 
of studies to Ce compounds with the AuCu3-type struc- 
ture: CePd3, CeRh 3 , Celn 3 , and CeSn 3 . Each of these 
materials is classified as a typical example of strongly 
correlated 4/ systems. 5, 14 ~ 19 ) 

CePd 3 is a typical heavy Fermion system with Tk of 
about 250 K, and has been studied extensively by various 
methods. 5 ' 18 ~ 31 ) It has a nonmagnetic Fermi liquid (FL) 
ground state at low temperatures. A controversy existed 
in the experimental works on inelastic magnetic excita- 
tion, but it has recently been resolved by a detailed study 
of the wave number vector (wave vector) dependence of 
spectra. 19, 21 ~ 23 ' The wave-vector-integrated-spectrum of 
the magnetic excitation has a broad peak at approxi- 
mately 55 meV. In the single-particle excitation, a broad 
peak with the 4/ character was observed on the inverse 
photoemission spectrum (IPES) side (i.e., in the energy 
region above the Fermi energy (-Ef)). This seems to be 
consistent with the high Tk of this compound. 24 ~ 27 ) How- 
ever, a strong peak structure has never been observed on 
the photoemission spectrum (PES) side (i.e., in the en- 
ergy region below Ef), contradicting other physical prop- 
erties. 28 ' 29 ' Recent careful study using the 3<i-4/ high- 
resolution resonant photoemission spectrum (RPES) re- 
vealed that the bulk component of the PES of this com- 
pound shows strong intensity at E-p consistently with 
high T K . 30) The angle-resolved PES (ARPES) has also 
been studied recently. 31 ' It may be worthwhile whether 
to confirm these recent results of studies are reproduced 
or not by the first-principles DMFT band calculation. 
CeRh 3 is known as one of the compounds having the 
most itinerant 4f states. 25, 26 > 32 ~ 36 ) 
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CeSri3 is also known to show the nonmagnetic 
FL ground state with a high characteristic tempera- 
ture. 37,38 ) Celn 3 has an antiferromagnetic ground state 
with a Neel temperature of Tn = 11 K. Recently, it was 
found that Tn decreases to zero under the pressure of 2.5 
GPa, in addition, the transition to superconductivity oc- 
curs at Tgc = 0.2 K. 39 ~ 41 ) It will be interesting to study 
the change of the band structure under pressures by the 
DMFT calculation. 

These compounds commonly have the AuCu3-type 
crystal structure, which is classified into the simple cu- 
bic (sc) lattice. The DMFT band structure of these com- 
pounds will be reported in two papers. In CePd3 and 
CeRh 3 , the hybridization of 4/ states with Ad states of 
the transition metal is very strong. The Ad states almost 
sink to below Ep in CePd3 whereas they are located near 
Ep in CeRh 3 . 42 ) In this paper, calculations for these Ad 
compounds will be reported. Calculations for CeSn 3 and 
Celn3 will be given in a subsequent paper. Their 4/ states 
hybridize with broad 5p states of ligand ions. 

Calculated results of CePd3 and CeRh 3 generally show 
reasonable agreement with experimental results of the 
PES, IPES, ARPES, and inelastic magnetic excitation 
by neutrons. However, the calculation gives a higher Tk 
than that expected from experiments when it is examined 
in detail. The DMFT band calculation for CePd 3 gives 
a Fermi surface (FS) structure similar to that obtained 
by the local density approximation (LDA) calculation at 
very low temperatures. When the temperature increases, 
4/ bands shift to the high-energy side and their lifetime 
broadening increases. This leads to the change of the FS 
structure into one that is similar to the FS of LaPd 3 . At 
T = 150 K, the FS has a different form from both FSs of 
LaPd3 and the LDA band of CePd3 as an intermediate 
stage of the change. The lifetime broadening overcomes 
the fine wave vector dependence of the 4/ spectrum at 
T = 300 K. In CeRli3, the 4/ band is located at about 
0.9 eV above Ep, and the dispersion of the DMFT band 
is almost identical to that of the LDA band in the en- 
ergy region near Ep. The density of states (DOS) has 
an appreciable value slightly below Ep in the ARPES. 
This low-binding energy part shows a weak wave vector 
dependence, though no flat bands do not exist in the 
vicinity of the Fermi energy. 

In §2, we briefly give the formulation on the basis of 
the LMTO method. Results of the application to CePd 3 
are shown in §3, and results for CeRh 3 are given in §4. A 
summary is given in §5. In the appendices, notes on the 
calculation of the total electron number are given. An 
efficient method of calculating the Cauchy integral using 
the spline interpolation scheme is also presented. This 
integral is frequently used in the DMFT calculation. 

2. Formulation 

The method of calculation is described briefly because 
its details have been given in previous papers. 6, 7 ) We 
consider the excitation spectrum of the following Hamil- 



tonian: 

n = n hDA + \ c ; . c ^ r7 - nf^f. (i) 

i r, 7 

Here, c^ a j r ^ is the annihilation operator for the atomic 
localized state 0| r (r ) at site i with the 7 orbital of 

the T-irreducible representation. The quantity n ^ A * is 
determined using the occupation number of the atomic 
4/ electron per Ce ion in the LDA calculation. We as- 
sume that the local Coulomb interaction acts only on the 
orbital ch a - „ . 

The excitation spectrum is expressed by introducing 
the self-energy terms, 6 ) 

"Hdmft = "Hlda 
+ (Sr( £ + i<5)+e r -e r DA )|^ r7 X^ r7 |, (2) 
i ,(r, 7 ) 

where ef; is the single-electron energy level of the 4/ 
state, and £p DA is the energy level in the LDA calcu- 
lation. The self-energy Sr(£ + ii5) is calculated by solving 
the auxiliary impurity problem with the use of NCA/ 2 vc; 
its outline is described in the Appendix of ref. 7. 

In later calculations we will approximate the localized 
4/ state <f) a by the band center orbital </>(—) because its 
localization is good for the 4/ state. </>(—) has the loga- 
rithmic derivative — I — 1 on the muffin-tin surface. 12,13 - 1 

In the LMTO method, the Hamiltonian "Hlda is diag- 
onalized using the LMTO bases, 

^H = E<a^). 0) 

qL 

Here, cr^ is the expansion coefficient of the jth eigen- 
vector on the LMTO base of the Bloch type, x^ L {r ), 
with the wave number vector k , the angular momentum 
(i,m), and the spin (a) at site q in the unit cell, where 
L = (f,m,a). 13 ) 

The explicit expression of X^,( r ) nas been given in a 
previous paper. 6 ) Note that they are not orthogonal to 
each other, but the eigenvectors ip^ are orthonormal. 

The DMFT band structure is calculated in the follow- 
ing way: (A) first the LDA part of the Hamiltonian 'Hlda 
is diagonalized for a given k , (B) then the matrix equa- 
tion of the Greenian is prepared in the manifold of the 
eigenvectors ip^ . The Greenian equation for the given 
k is written as 

[zI-D Tj - DA (k)-X(z)]G(z;k)=I, (A) 

where / is the unit matrix and -Dlda (k ) is the diagonal 
matrix of the eigenenergies of Hlda with k . The matrix 
elements of S(z) are given by calculating the self-energy 
operator term of eq. (2) based on cq. (3). 
The DOS on the atomic 4/ state is given by 

p r hand \e; k ) = --^tr[d r G(e + id; k )], (5) 
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where the projection operator is defined as 
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The local DOS in the DMFT band calculation is obtained 

over k in the Brillouin zone: 

37 Efc Pr and) ( e ; fe )- 43) Here ' ^ is the total 
number of unit cells. 

The auxiliary impurity problem is solved by the 
NCA/ 2 vc method. The splitting of the self-energy due 
to the SOI and CFS effects is considered. As shown in 
ref. 7, this method gives an accurate order of the Kondo 
temperature when the result is compared with that of the 
more correct numerical renormalization group (NRG) 
calculation 7 ' 44 - 1 in a simple model case. 

Since the method of the self-consistent calculation in 
the DMFT has been described previously, 6, 7 > we ex- 
clude the detailed explanation from this paper. First of 
all, we calculate the self-consistent LDA band by the 
LMTO method, and potential parameters, except for the 
/ levels, are fixed to those in the LDA calculation. (I) 
We calculate the atomic 4/ density of states p r imp ' ) (e) 
(4/ DOS) for the auxiliary impurity Anderson model by 
the NCA/ 2 vc method with a trial energy dependence of 
the hybridization intensity (HI) and 4/ levels, 45 > then 
calculate the local self-energy. (II) The DMFT band cal- 
culation is carried out using the self-energy term, and 
the local 4/ DOS in the DMFT band is calculated. The 
calculation is iterated so that the 4/ DOS of the local 
auxiliary impurity model and the DMFT band satisfy 
the self-consistent conditions. 46, 47 * ) 

The 4/ level is adjusted in the DMFT self-consistent 
iterations under the condition that the 4/ occupation 
number has a given target value, (rsl. target), which is 
estimated from the LDA band calculation. The tempera- 
ture dependence of the Fermi energy, Ep , is neglected by 
fixing it at a value determined at a low temperature. It is 
estimated using the occupation number of the renormal- 
ized band (RNB) calculation, in which the self-energy is 
approximated by an expansion form up to the linear term 
in the energy variable at Ep (see Appendix A for the cal- 
culation of the total occupation number). The target 4/ 
electron number rif (rsl. target) is imposed on the occu- 
pation number calculated directly using the resolvents 
to stabilize the self-consistency iterations. The occupa- 
tion number obtained by the integration of the 4/ DOS, 
rif (intg.), has a deviation within 1.0% from rif (rsl. target) 
because many intermediate calculation process are in- 
cluded. 

3. CePd 3 

3. 1 Density of states 

In Fig. 1, we show the DOS of the 4/ component 
(4/ DOS) for CePd 3 at T = 37.5 K. The solid line shows 
the total 4/ component of the PES (4/ PES). The dashed 
line is the DOS of the T-j component and the dot-dash 
line is the DOS of the Tg component. The two-dots-dash 
line is the DOS of the j = 7/2 component. The CFS of 
the self-energy in the excited j = 7/2 multiplet is ne- 
glected. The vertical dot-dash line indicates the Fermi 




0.10 0.30 0.50 0.70 0.90 1.10 1.30 
energy (Ry) 

Fig. 1. 4/ DOS of CePd 3 at T = 37.5 K. The solid line shows the 
total 4/ PES. The dashed line is the DOS of the (j = 5/2)F 7 
component, the dot-dash line is the DOS of the (J = 5/2)rg 
component, and the two-dots-dash line is the DOS of the j = 7/2 
component. The Fermi energy E-p = 0.6143 Ry is indicated by 
the vertical dot-dash line. The inset shows spectra in the vicinity 
of Ep. 
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0.10 0.30 0.50 0.70 0.90 1.10 1.30 
energy (Ry) 

Fig. 2. Hybridization intensity (HI) 54 ) of CePd 3 at T = 37.5 K. 
The dashed line is the HI of the (j = 5/2)T7 component, the 
dot-dash line is the HI of the (j = 5/2)rg component, and the 
two-dots-dash line is the HI of the j = 7/2 component. Ep is 
indicated by the vertical dot-dash line. The inset shows HI in 
the vicinity of Ep . 



energy E-p = 0.6143 Ry. The inset shows spectra in the 
energy region near Ep . The 4/ DOS has a large peak at 
0.642 Ry, about 0.028 Ry (0.38 eV) above Ep. This has 
mainly the j — 7/2 character. The spin-orbit splitting 
on the IPES side is usually enhanced in Ce systems. The 
spectral intensity on the PES side consists mostly of the 
component. The PES has a sharp peak at Ep with a 
steep tail up to 0.094 Ry (1.3 eV) below Ep. It also has 
a long tail with small structures reflecting the DOS of 
Ad states of Pd on the high-binding-energy side. In the 
steep tail region, shoulders appear at binding energies of 
0.004 Ry (0.05 eV) and 0.024 Ry (0.33 eV). These may 
correspond, respectively, to the CFS and the SOI side 
band. The sharp peak at Ep with the steep and the long 
tail has been observed in high-resolution experiments by 
Kasai et al. 30 ^ The shoulder due to the SOI seems to be 
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0.10 0.30 0.50 0.70 0.90 1.10 1.30 
energy(Ry) 

Fig. 3. Total DOS (dashed line) and 4/ DOS (solid line) of CePd 3 
at T = 37.5 K. The spectra are broadened by an imaginary 
factor, 7 = 0.01 Ry, in the energy variable. Ep is indicated by 
the vertical dot-dash line. 



Table I. Various quantities obtained in the DMFT calculation for 
CePd3 at T = 18.75 K. n[! mp is the occupation number in the 
auxiliary impurity problem of the effective HI of the DMFT cal- 
culation, er is the energy level in Ry. pp(Ep) is the partial DOS 
at Ep in Ry - 1 . 1 = 1 — d^FEp (s) / ds\ e f is the mass renor- 
malization factor of the 4/ band, ep is the effective energy of 
the renormalized band in Ry, and Tp is the imaginary part of 
the self-energy in Ry. The effective energy levels are measured 
from the Fermi energy Ep = 0.6143 Ry. The lattice constant 
is a = 7.80079 a.u., and the spin-orbit interaction constant is 
Cif = 7.2302 X 10~ 3 Ry. The 4/ level in the band calculation 



'4/ 



0.6754 Ry. The electrostatic CFS of {j = 5/2)r 7 is 



set to be AE r7 (= ef 



= 150 K. The target 4/ occu- 
pation number is n t (rsl. target) = 0.98, and the resultant oc- 
cupation number calculated using the resolvent is 0.980. The 
4/ occupation number calculated by integrating the spectra is 
rij>(intg.) = 0.989, and the obtained total band electron num- 
ber is 7V(total; RNB) = 34.006. E mc \ is the characteristic en- 
ergy of the quasi-elastic excitation and i?CFS is the CFS exci- 
tation energy estimated from the magnetic excitation spectrum. 
The occupation number of the 4/ state in CePd3 in LDA is 
rif(Ce, LDA) = 1.045. The Coulomb constant U is set to be 
0.51 Ry (6.9 eV). 45 ) 





r 7 


r 8 


3 = 7/2 


n r imp ' ) 


0.084 


0.719 


0.185 


£r(Ry) 


-0.19027 


-0.19182 


-0.16493 


PriEp^Ry- 1 ) 


3.8 


59.0 


0.5 




6.3 


6.7 


3.6 


£r(Ry) 


0.0153 


0.0057 


0.1184 


Fr(Ry) 


2.7 xl0~ 4 


16.0 xlO" 4 


3.0 xl0~ 4 


B ine i = 27 meV, 


EcFS = 


41 meV 





identified. 

The Kondo temperature 7k and the CFS excita- 
tion energy are, respectively, estimated to be about 
27 meV (310 K) and 41 meV (480 K) from magnetic 
excitation spectra, as will be shown in Fig. 4 in the next 
subsection. We note that spectra do not show apprecia- 
ble change even when the temperature is decreased to 
18.75 K in the calculation. 

The parameters and the calculated values are given in 
Table I. The LMTO band parameters 13 ^ for states ex- 
cept for the / component are fixed to those of the LDA 



calculation. Ep is fixed to the value determined by the 
occupied state in the RNB, as is discussed in a later 
section. The relative occupation number of the r 7 com- 
ponent to the r 8 , 0.084/0.719 (0.12) is small compared 
with the 0.5 expected from the ratio of the degeneracy, 
but is very large compared with the value expected from 
the simple model of the CFS for an isolated ion with 
e r7 - er8 = 0.00155 Ry (240 K). The ratio does not 
change so greatly even when the temperature is raised: 
0.089/0.723 (0.12) at 150 K and 0.108/0.707 (0.15) at 
300 K. Moreover, the occupation of the j = 7/2 compo- 
nent, 0.185, is not small, and is almost independent of T. 
A simple picture of the CFS for an isolated ion cannot 
be applied. 

Usually, the electrostatic potential causes cubic CFS 
in 4/ electron systems with a higher energy level of 
about 150 K for the T 7 state. 48 ' 49 ) This is included in 
the present calculation. Even when it is neglected, the 
4/ DOS and the magnetic excitation spectrum are not 
greatly changed because the hybridization effect causes 
large CFS, of greater than 300 K. 50 ) 

The target value of the occupation number on the 
atomic 4/ states, rif (rsl. target), has been tentatively cho- 
sen to be 0.98 in the present calculation. This value is 
small compared with the 4/ occupation number, 1.045, 
of the LDA band calculation for CePd3. When we carry 
out the LDA calculation for compounds in which Ce 
ions are replaced by La ions, the occupation number 
on the 4/(La) state usually amounts to O.l. 51 ' We use 
the occupation number as reduced to 94% of the LDA 
value as the atomic occupation number on the Ce 4/ 
state. Kanai et al. concluded that 4/ occupancy in 
CePd3 is expected to be 0.92 on the basis of the re- 
sults of resonant inverse photoemission (RIPES) experi- 
ments. 27, 52 > 53 ) When we perform the DMFT band calcu- 
lation by setting n/ (rsl. target) to be 0.92, 7k determined 
from the magnetic excitation spectrum is expected to be 
500 K. On the other hand, 7k is greatly reduced to about 
10 K when we choose 1.05. 

In Fig. 2, we show the trial HI 54 ) obtained in the 
DMFT band calculation at 37.5 K. It has very large 
peaks corresponding to the Ad band of Pd, but these 
peaks are located in the energy region deep below Ep. 
The HI is not high in the energy region near Ep . In par- 
ticular, the HI of the r 7 component is low though it is 
high in the deeper energy region. The overall features of 
the HI in the DMFT band are similar to those calcu- 
lated directly using the LDA band, but the HI in DMFT 
is increased in the vicinity of Ep to about twice the LDA 
value for the T 7 and Tg components, and is decreased 
in the deep energy region. Moreover, the DMFT calcu- 
lation causes fine structures of the HI in the vicinity of 
Ep, which are shown in the inset of the figure. This con- 
trasts with the HI in the LDA, which has a weak energy 
dependence in this region. The HI of the Tg component 
in DMFT has a small peak at Ep, while those of the T 7 
and j = 7/2 components have small peaks above or on 
both sides of Ep. The reason for this different behavior 
is not clear at present, but we note that the peaks of 
the 4/ DOS for the latter two cases are located above 
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Ep. When we do a calculation in the fictitious case that 
Tj is mainly occupied by assigning a low energy level to 
it, the HI of T7 has a small peak at Ep. However, this 
result should not be used as a general rule, because the 
modification of HI in DMFT is delicately dependent on 
details of the band structures. The Kondo temperature 
is increased in the DMFT calculation in the CePd3 case. 
We have obtained the Kondo temperature of 10 K in 
the single impurity calculation using the HI of the LDA 
band. 

The total DOS at T = 37.5 K is shown by the dashed 
line in Fig. 3. The large peaks at about 0.3 and 0.45 
Ry have the 4d character of Pd, and that at 0.8 Ry has 
the 5d character of Ce. These peaks are also obtained 
by the LDA calculation. 55 ' The sharp 4/ peaks slightly 
above Ep, which are called the f 1 peak in IPES, are 
also obtained in the LDA calculation. 42 ' Their intensity 
in DMFT is reduced compared with that in the LDA 
because a part of it is transferred to the intensity of the 
peak at about 1.0 Ry, which is called the f 2 peak. In 
the present calculation, the width of the f 2 peak is not 
large because the multiplet splitting of the f 2 final state 
is neglected. 

In the analysis of RIPES experiments, the ratio of the 
f 1 peak to the total RIPES intensity has been given 
as 0.22, 27 ) whereas it is estimated to be about 0.2 in 
the present calculation. The present DMFT calculation 
seems to give results that emphasizes the hybridization 
effects strongly (i. e. the higher 7k). 

3.2 Magnetic excitation 

-120meV 120meV 

— I 1 h- 

300 i I 




-0.01 0.00 0.01 

E (Ry) 

Fig. 4. Magnetic excitation spectrum of CePd3. The solid line 
shows the spectrum at T = 37.5 K, the dot-dash line is the spec- 
trum at T = 150 K, and the two-dots-dash line is the spectrum at 
T = 300 K. The dashed line is the spectrum in a fictitious case 
where matrix elements of the magnetic moment are restricted 
within the intra-Tg manifold of space. 

In Fig. 4 we show the magnetic excitation spectrum. 
The total magnetic excitation spectrum at T = 37.5 K is 
shown by the solid line. It has a peak at about E = 0.003 
Ry (41 meV) . The dashed line depicts the spectrum for a 



fictitious case in which the matrix elements of the mag- 
netic moment are nonzero only in the manifold of Tg, and 
thus it may correspond to the excitation spectrum within 
the r 8 manifold. It has a peak at about E = 0.002 Ry (27 
meV), and leads to a faint shoulder in the solid line. The 
CFS excitation energy seems to be slightly larger than 
7k in the present calculation. The calculated magnetic 
susceptibility is 1.9 xlO -3 emu/mol, whereas the exper- 
imental value is 1.5 xl0~ 3 emu/mol at T ~ K. 56 ' 57 ' 
We show magnetic excitation spectra at T = 150 K 
and at T = 300 K by the dot-dash line and the two- 
dots-dash line, respectively. The spectrum at T = 300 K 
has a broad peak centered at about 0.003 Ry (41 meV). 
The shoulder shifts to E ~ and becomes a peak, the 
overall features do not change so greatly when we ne- 
glect the electrostatic CFS of 150 K. The magnetic exci- 
tation spectrum observed in the wave-vector-integrated 
case has a peak at about 55 meV at T = 10 K, and the 
peak shifts to the low-energy side as T increases. 21 ' The 
calculated results seem to be generally consistent with 
those of the experiment. However, the peak at E ~ is 
higher at T = 300 K in the experiments. 21 - 1 The present 
DMFT calculation seems to give stronger HI compared 
with the value in the experiments. The detailed calcula- 
tion of physical quantities using finely tuned parameters 
will be given in the future. 

In ref. 22, excitation peaks with an energy of 15 meV 
and less than 3 meV are indicated. Low-energy peaks 
are not expected in the present calculation of the wave- 
vector-integrated spectra because 7k is not low. One pos- 
sibility of the origin of the low-energy peaks may be the 
wave vector dependence of the magnetic excitation spec- 
tra, as noted in ref. 23. 

3.3 RNB, and wave number-vector-dependent DOS 

In Fig. 5, we show the RNB dispersion at T = 37.5 
K. The energy shift (the real part of the self-energy at 
Ep: 5RSr(£ ; F)) and the mass renormalization factor (1 — 
d^fEr(e)/de\E F ), which are given in Table I, are taken 
into account in this calculation. Narrow bands with the 
j = 5/2 character appear slightly above Ep, and those 
with the character of j = 7/2 appear around the energy 
0.645 Ry which is near the energy of the j — 7/2 peak 
in the 4/ DOS shown in Fig. 1. 

The lowest 4/ band sinks to below Ep near the T and R 
points, and is located above Ep in the other regions. The 
dispersion of RNB corresponds well with the behavior of 
the fc -dependent density of states (k-DOS). For example, 
we show the k-DOS when k moves from the T (bottom) 
to the R (top) point along the A line in Fig. 6. A peak of 
the DOS with mainly the 4/ character is located below 
Ep at the T and R points. Starting from the peak below 
Ep at the T point, one of ridge lines runs above the Fermi 
energy across Ep, and then connects to the peak below 
Ep at the R point. Another runs to the low-energy side of 
Ep up to 0.611 Ry at the halfway, and then turns back 
to connect to the peak below Ep at the R point. This 
"hanging" branch does not cross the Fermi energy. 

Note that we have depicted the total spectral inten- 
sity, not the /-component, in Fig. 6. In the energy region 



6 J. Phys. Soc. Jpn. 



Full Paper 



Osamu Sakai 



20 




Fig. 6. Wave number vector (k ) dependence of the DOS (k-DOS) 



for CePd3 at T = 37.5 K. k is moved from the T point (bot- 




Fig. 5. Band dispersion of the renormalized band (RNB) picture 
for CePd3 at T = 37.5 K. The symbols under the horizontal axis 
denote the symmetry points and axes of the Brillouin zone of the 
sc lattice. E-p is indicated by the horizontal dashed line. 






0.61 0.62 
energy (Ry) 

Fig. 7. k-DOS of CePd 3 at T = 150 K. k is changed from the T 
point (bottom) to the R point (top) along the A axis. 



shown in the figure, the spectral intensity has mainly 
the 4/ character. On the other hand, very sharp spectral 
peaks of a non-4/ character appears in the energy region 
outside of the figure. The hanging band on the A axis has 
a stronger non-4/ character. This branch is a hybridiza- 
tion band between the 4/ of Ce and a conduction band 
that has the character of the sp-free electron band and 
the 5d of Ce. 

The dispersion of the RNB is qualitatively similar to 
the result of the band structure determined by Hasegawa 
and Yanase by the LDA calculation, 5S ) although the 
width of the 4/ band with j = 5/2 is about 0.005 Ry 
in the RNB, whereas that of the LDA band is about 
0.02 Ry. Both calculations give electron pockets at the T 
and R points, and hole pockets centered on the T axis. 59 ) 

The Main features of the band dispersion near the 
Fermi energy are formed by the hybridization of the nar- 
row 4/ bands and the wide sp bands. The 4/ bands have 
dispersion characterized by the LMTO (linear combina- 
tion of atomic orbitals) tight-binding bands of the sc lat- 
tice. Although the 4/ band width is reduced in DMFT, 
the qualitative features of the dispersion of the 4/-sp hy- 
bridized bands are not changed because the number of 



participating sp bands is small and their dispersion is 
very rapid compared with that of 4/ bands. 

Here, we must note a weak point of the NCA/ 2 vc 
method, that is it does not automatically ensure the 
Fermi liquid sum rule, i.e., the integral of the total 
DOS below E F , iV(total; p) = J de^J^k (-^ G ( £ + 
id; k ))/(e), is not equal to the occupation number of elec- 
trons calculated by the volume of the occupied states in 
the k space. The quantity N (total; p) is a smaller value 
in CePd3, and thus we obtain a higher Fermi energy. If we 
use it in the RNB calculation, the volume of the occupied 
states in the k space becomes large. In the case of CePd3, 
which has an even total electron number, the balance be- 
tween the electron and hole states is lost. In this study, 
we tentatively use the Fermi energy determined using the 
quantity N (total; RNB) = i J2\k /(^RNb(AA; )) to en- 
sure the electron-hole balance in the RNB band, where 
£ , RNB(Afe ) is the energy obtained by the RNB calcula- 
tion. (For more details, see Appendix B.) 

When we calculate the RNB dispersion at T = 150 K, 
the 4/ band shifts slightly to the high-energy side. The 
4/ state at the T point nears Ep, and the 4/ state at 
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Fig. 8. Angle-resolved PES of CePd 3 at T = 37.5 K. The surface 
is assumed to be (111). The representative fc is swept from the 
r (bottom) point to the M (top) point along the S axis, and 
the intensities are averaged over the wave number vector normal 
to the (111) surface. This means that fc also runs from the R 
(bottom) point to the X (top) point along the S axis. The Fermi 
energy E-p is indicated by the vertical dot-dash line. The total 
DOS is shown, but the intensity above the energy of 0.612 Ry 
is dominated by the 4/ character. The intensity below 0.611 Ry 
has mainly the non-4/ character. 



10 
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Fig. 9. Angle-resolved PES of CePd 3 at T = 150 K. For the def- 
inition of lines, see the caption of Fig. 8. 

the R point is located on Ep. Therefore, a part of the 
hanging band on the A line rises above Ep. Hole pockets 
on the T axis grow into larger hole regions around the R 
point. As the temperature increases further, the 4/ state 
at the R point shifts up to above Ep, and a large hole 
surface enclosing the R point appears. In other words, 
we have a connected electron Fermi surface (FS) that 
contains the X and M points inside it. 

The primary structure of the FS of the RNB at high 
temperatures is similar to that of the FS of the LDA 
band of LaPd3, but their fine topologies will differ. At 
T = 150 K, the main part of the hanging band on the A 
line is still located below E-p. The hole sheets around the 
r and R points are separated by an electron region on 
the A axis. On the other hand, this entire hanging branch 
is located above Ep in LaPd 3 , and thus the separation 
by the electron region does not occur. 

We should, of course, note that the RNB picture has 



only limited meaning at high temperatures because the 
imaginary part is large. 60) The k-DOS at T = 150 K for 
fc along the A line is shown in Fig. 7. The widths of peaks 
become large, but we can see that the peak at the R point 
is located almost at Ep, and the peak at T nears Ep. The 
trace of broad peaks shows a shift corresponding to that 
of the RNB dispersion. We can recognize that the ridge of 
the peaks of the hanging band crosses the Fermi energy, 
as noted in the previous paragraph. At T = 300 K, the 
width of peaks becomes so large that it surpasses the fine 
fc dependence of the spectra. However, the 4/ peak near 
the Fermi energy still exists in the 4/ DOS calculation, 
similar to that shown in Figs. 1 and 3. It moves slightly 
to the high-energy side with increasing width at 300 K. 

The present calculation has been performed with fixed 
Ep and n/(rsl. target). The calculated total occupation 
number in the DMFT band increases to N {total; RNB) = 
34.304 at T = 300 K, although the energies of the 4/ 
bands shift upward. 61 ) We point out that even when we 
move E-p, the energy of 4/ bands relative to the Fermi 
energy will not change greatly because the Kondo res- 
onance peak usually shifts following the change of the 
Fermi energy. 

Recently, the wave vector dependence of the PES has 
been observed under the 4d — >• 4/ resonant condition. 31 ) 
The experiment was carried out for the (111) surface of 
a thin film by sweeping the fc vector along T-K-M-K-r 
in the surface Brillouin zone (BZ). This may correspond 
to the following sweeping of fc in the 3-dimensional sc 
BZ: the component parallel to the surface runs as r-(E)- 
M-(S)~r with averaging over its components normal to 
the surface. This means that the sweeping of the parallel 
component R-(S)-X-(S)-R is also included. 

In Fig. 8, we show the PES when the representative 
wave vector moves from the T (bottom) point to the M 
(top) point with average intensities for fc normal to the 
surface. The total intensity is plotted in the figure, but 
the intensity above the energy of 0.612 Ry has mainly 
the 4/ character. The spectra have a peak for fc near 
r, and this corresponds well to the experimental results. 
The intensity below the energy of 0.611 Ry has a non4/ 
character. As noted previously, the contribution from the 
R point is also included at the representative T point. At 
these points, the 4/ component is located below Ep at 
T = 37.5 K. In Fig. 9, we show spectra at T = 150 
K. The peak of the intensity at T nears Ep because of 
the shifting up of the 4/ band. The thermal distribution 
effect and the very high intensity of the 4/ DOS above 
E-p also have affect this spectral shape. The change of 
the peak can be checked in experiments. 

4. CeRh 3 

In Fig. 10, we show the 4/ DOS of CeRh 3 at 10 3 K. 
The 4/ spectrum has a sharp and large peak at an en- 
ergy about of 0.76 Ry, which is 0.066 Ry (0.9 eV) above 
the Fermi energy. This separation of energy from Ep is 
slightly smaller than that of the 4/ band in the LDA, 
about 0.009 Ry (1.2 eV). The result seems to be consis- 
tent with those of RIPES experiments and their detailed 
analysis. 26 ' 32 - 1 The PES has a relatively large peak at 
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Fig. 10. 4/ DOS of CeRh 3 at T = 10 3 K. For the definition of 
lines, see the caption of Fig. 1. E-p = 0.6940 Ry is indicated by 
the vertical dot-dash line. 
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Fig. 12. Angle-resolved PES of CeRh 3 at T = 300 K. For the 
definition of lines, see the caption of Fig. 8. The k point corre- 
sponding to the second solid line from the top is named fc * in 
the text. 
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Fig. 11. Total DOS (dashed line) and 4/ DOS (solid line) of 
CeRh 3 at T = 10 3 K. See the caption of Fig. 3. 

about 0.1 Ry (1.4 eV) below Ep and a peak at E-p. The 
qualitative behavior of the present PES is similar to the 
result obtained by Harima in the LDA. 62 ) However, the 
binding energy of the peak at 0.1 Ry below Ep is slightly 
lower and the intensity is higher than that of the LDA. 
The peak at Ep is also slightly sharper. In the experi- 
ment, the sharp peak at Ep was observed, but the peak at 
0.1 Ry below Ep has not been identified at present. 33 ~ 36 ) 
In Fig. 11, we show the total and 4/ DOS of CeRh 3 
at T = 10 3 K. The Ad band of Rh exhibits strong peaks 
of the DOS at about 0.4 and 0.6 Ry The 4/ component 
also has a small peak at about 0.6 Ry, as noted previ- 
ously. The Fermi energy is located in the top region of 
the Ad band, and the hybridization intensity in this re- 
gion is high. Dispersions of the RNB are almost identical 
to those of the LDA in the energy region very near Ep, 
but the width of the 4/ band, which is located at 0.76 R, 
is about 2/3 that of the LDA. A f 2 satellite peak ap- 
pears in DMFT at about 1.2 Ry on the high-energy side. 
The ratio of the intensity of peak at 0.76 Ry (f 1 peak) 
to the total IPES intensity is estimated to be about 0.4 
in the present calculation, while a slightly larger value, 
0.6, was obtained in the experiment. Uozumi et al. pre- 



Fig. 13. Magnetic excitation spectra of CeRh 3 at T = 10 3 K. 
The solid line is the total intensity and the dashed line is the 
result calculated using a fictitious model in which the matrix 
elements of the magnetic moment are restricted within the j = 
5/2 manifold of states. 

dieted the 4/ occupation number to be 0.86, 32 * 1 but we 
tentatively used n/(rsl. target) = 0.94, which is 94% of 
the LDA value. The difference between these values is 
not small, but will not cause extreme differences in the 
calculation of the very strong hybridization limit. The 
mass enhancement factor is expected to be about 2, as 
given in Table II. 

In Fig. 12, we show the PES for the (111) surface when 
the representative wave vector moves from the Y point 
(bottom) to the M point (top). In the energy region be- 
tween Ep and 0.6 Ry, the relative intensity of the 4/ 
component is about 10% of the total. The spectra have 
fine peaks as if some flat bands exist slightly below Ep. 
However, we note that no flat bands exist very near Ep 
in the RNB dispersion. Let us denote as k * the represen- 
tative k corresponding to the second solid line from the 
top (i.e., the k point on the E axis at a distance of about 
0.3 x length of the E axis from the M point). The spectral 
intensity slightly below Ep is relatively large for k near 
k * . For the wave vectors near k * , several bands stay 
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Table II. Various quantities obtained in the DMFT calculation 
for CeRli3 at T = 10 3 K. For the definition of notation, see the 
caption of Table I. The Fermi energy is Ep = 0.6940 Ry and 
AE7 = K. The lattice constant is a = 7.6005 a.u., and the 
spin-orbit interaction constant is £4^ = 7.3374 X 10~ 3 Ry. The 
4/ level in the band calculation is ^4p A = 0.7639 Ry. The target 
4/ occupation number is n^(rsl. target) = 0.94, and the resultant 
occupation number calculated using the resolvent is 0.940. The 
4/ electron number calculated by the integrating the spectra is 
n^(intg.) = 0.943, and the obtained total band electron number 
is A^total; RNB) = 30.994. n^(Ce, LDA) = 0.996. 
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slightly below Ep around the M and X points when the 
normal components are varied. Calculated results show 
quantitatively similar behaviors to experimental results, 
but the careful separation of the surface and bulk com- 
ponents is necessary to enable a detailed comparison. 63 ) 

The magnetic excitation spectra are shown in Fig. 13. 
They have a steep increase at the excitation energy of 
about 0.03 Ry. This energy may correspond to the en- 
ergy from Ep to the low-energy edge of the peak at 0.76 
Ry of the 4/ DOS. The-low energy end of the spectrum 
mainly originates from the excitation within the j = 5/2 
components, but contribution of the j — 7/2 components 
is not small even in the low-excitation-energy region. The 
j = 7/2 components also join the Kondo effect. We may 
expect the Kondo temperature of this system to be about 
0.03 Ry (4700 K). The calculated value of the magnetic 
susceptibility is 0.7 x 10~ 4 cmu/mol. The experimental 
magnetic susceptibility (x~4x 10^ 4 emu/mol) 57 ' may 
indicate a high Tk of several thousand K, but the calcu- 
lated value of 4700 K seems to be too high. DOSs, both 
of the total and of the 4/ component, do not have any 
gap, as seen in Fig. 10, but the magnetic excitation spec- 
trum has a shape indicating the existence of a pseudogap. 
We note that, in Fig. 13, the excitation spectra of only 
intra-4/ components are shown. A broad continuous 5d 
component will be superposed on these spectra. 

The HI of CeRli3 in the LDA calculation has a spec- 
trum shape similar to that of the partial DOS on the 5d 
state of Rh, i.e., the spectral shape given by subtracting 
the 4/ and 5d parts of the Ce ion from the total DOS in 
Fig. 11. The HI in DMFT is almost equal to that of the 
LDA in the high-energy region, but it has a steep dip at 
an energy slightly above Ep. Similar behavior has been 
seen in the HI of the j = 7/2 component of CePd3 which 
has shown in Fig. 2. 

5. Summary and Discussion 

We have studied the electronic structures of CePd3 and 
CeRli3 on the basis of the DMFT calculation. The aux- 
iliary impurity problem was solved by a method named 
NCA/ 2 vc, which includes the correct exchange process 



of the f 1 -> f° and f 1 — > f 2 virtual excitation. The 
splitting of the self-energy owing to the SOI and CFS 
effects was also considered. 

The DMFT band calculation gives Fermi surface struc- 
tures similar to those obtained by the LDA calculation 
in CePd3 58 ) at very low temperatures. Electron pockets 
appear at the bottom of the 4/ band at the T and R 
points. Hole pockets appear that are centered on the T 
symmetry axis. 

The 4/ band shifts to the high-energy side relative to 
the Fermi energy as the temperature increases. At the 
same time, the lifetime width of the 4/ states increases. 
The band structures produced by the band overlap be- 
tween the 4/ and non-4/ components shift up to the high 
energy side of the Fermi energy. Therefore the primary 
structures of the band in the vicinity of the Fermi energy 
approaches to those of the sp free-electron-like band of 
LaPd 3 . However, some characteristic features of the LDA 
band of CePd3 remain at higher temperatures. For ex- 
ample, a region that electrons occupy will appear on the 
A axis in CePd3 at T ~ 150 K, which is not expected in 
LaPd3. In CcPd3, the lifetime broadening overcomes the 
k dependence of the 4/ spectrum at room temperature, 
thus, the 4/ state becomes a broad dispersionless state 
located above Ep. 

The ARPES of CePd3 shows strong intensity near the 
representative T point at low temperatures because the 
4/ band is located below the Fermi energy at the T and R 
points. This result seems to be consistent with the recent 
experimental results. 31 ) The observed 4/ component will 
be greatly reduced at room temperature because of the 
shift of the 4/ band to the high-energy side. 

The k-integrated magnetic excitation spectrum has a 
peak at 41 meV and a faint shoulder structure at about 
27 meV at low temperatures. The temperature depen- 
dence of the excitation spectrum generally seems to be 
consistent with the results of experiments. 2 ^ The mag- 
netic CFS excitation energy is estimated to be about 
41 meV, while the Kondo temperature is slightly smaller, 
320 K (27 meV). 

The calculated PES shows good correspondence with 
the bulk component obtained in recent high-resolution 
experiments. 30 ^ The intensity ratio of the f 1 peak of the 
IPES to the total IPES is estimated to be about 0.2, 
while it was predicted to be 0.22 by the recent experi- 
ment analysis. 27 ^ The HI is enhanced near the Fermi en- 
ergy in the DMFT band compared with that of the LDA 
calculation. The present DMFT calculation seems to give 
electronic structures with slightly stronger HI than that 
expected from experiments for CePd3. 

The DMFT calculation for CeRh 3 gives an almost 
identical band structure to that obtained by the LDA 
calculation. However, the energy of the f 1 peak of the 
IPES in the former is slightly lower than that in the 
latter. In addition, the f 2 satellite peak appears in the 
DMFT band calculation. The calculated intensity ratio 
of the f 1 peak to the total IPES, 0.4, is comparable to, 
but smaller than the experimental value of O.6. 32 -* The 
calculated PES has a sharp peak at the Fermi energy. 
A peak reflecting the 4d(Rh) DOS also appears, simi- 
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larly to the result of the LDA calculation. The former 
has been observed, but the latter has not been identified 
in experiments. 35, 36 - ) 

The Fermi energy is located in the energy region of the 
Ad band of Rh in CeRh 3 , therefore the HI near the Fermi 
energy is strong, about three times greater than that of 
CePd3. The 4/ band width of CeRli3 is about twice that 
of CePd 3 in the LDA calculation. In the DMFT calcula- 
tion, the characteristic energy scales are drastically dif- 
ferent from each other, 300 and 5000 K. 

In CcRh 3 , the ARPES shows an appreciable DOS 
slightly density below the Fermi energy and is weakly 
dependent on the wave vector, although no flat 4/ bands 
exist near the Fermi energy in the RNB dispersion. The 
spectra have stronger intensity halfway along the £ axis 
from the T to M points of the Brillouin zone. This seems 
to be similar to the results of experiments, but further 
studies to separate the surface effects are necessary. 63 ) 

General features of the experimental results for CePd 3 
and CeRh 3 are reproduced by the DMFT band calcu- 
lation with the LMTO+NCA/ 2 vc scheme. However, the 
present DMFT calculation gives a higher Kondo tem- 
perature than that in the results of the detailed anal- 
ysis of experiments. Moreover, there is some arbitrari- 
ness in the choice of the target value of the 4/ occupa- 
tion number, n/(rsl. target). According to calculations for 
AuCu 3 -typc Ce compounds, accurate calculated results 
seem to be obtained when a 4/ occupation number be- 
tween 90% and 95% of the LDA value is used. We have 
used 94% in the present calculation, and we obtained 
310 K (27 mcV) for T K of CePd 3 (n/(rsl.target) = 0.98) 
and 4700 K (0.03 Ry) for T K of CeRh 3 (n/(rsl.target) = 
0.94) from the magnetic excitation spectrum. When we 
perform calculations using 90%, the 7k are 400 K for 
CcPd 3 (n/(rsl.target) = 0.94) and 6300 K for CeRh 3 
(n/(rsl. target) = 0.90). These values are not greatly dif- 
ferent from the previous values because these compounds 
belong to a group of materials having high 7k . For mate- 
rials with lower 7k, 7k drastically depends on the choice 
of n f (rsl. target) . Careful treatment of the target value is 
necessary in such cases. 

The effectiveness and some of the weaknesses (for ex- 
ample, the correct calculation of the occupied electron 
number) of the present DMFT scheme are recognized. 
Calculations of the 4/ band state in Celn 3 and CeSn 3 , 
and also various Ce compounds will be carried out in the 
near future. 

At the end of this paper, we refer the very early and 
recent application of methods similar to NCA/ 2 vc to 
the DMFT band calculation. Lasgsgaard and Svane cal- 
culated the band structure of Ce pnictides in 1998. 64 ^ 
Recently, Haule et al. studied the a — > 7 transition of 
Ce metal, 65 ) and Shim et al. studied the electronic band 
structure of Celrln 5 . 66 ) The CFS of the self-energy was 
not considered in those studies. 
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Appendix A: Occupation number 

The calculation based on the NCA/ 2 vc method usu- 
ally does not satisfy the FL sum rule, i.e., the integral 
of the total DOS below the Fermi energy, iV(total; p) = 

I dsjjJ2k ( — i^^( e + ))/( e )' i s not cc l ua l to the 
value N(total-lnG) = / de^ J^k (-^Slndet G{e + 

i<5; k where /(e) is the Fermi distribu- 

tion function at T = 0. In the RNB calculation, 
in which the self-energy term is approximated by 
E r (s) » 5RS r (S F ) + ™§^1\ Ef ( z - E F ), we obtain 
real eigen energies £?RNB(Afc ) where A is the band suf- 
fix. The number of occupied states in the RNB band, 
N (total; RNB) = i J2k /(^RNb(AA; )), is expected to 
agree with AT (total; In G) at T = if the imaginary part 
can be neglected near the Fermi energy. However, the 
imaginary component has a considerable magnitude in 
NCA/ 2 vc even at very low temperatures; this is partly 
because T must maintain the condition T > 0.17k- For 
example, AT (total; p), N (total; In G), and N (total; RNB) 
are 33.59, 34.16, and 34.01, respectively, at T = 18.75 K 
in CePd3. The difference between these values becomes 
serious when the detailed structure of the Fermi surface is 
discussed. In this study, we use ./V (total; RNB) to deter- 
mine the Fermi energy Ep at T = because this quantity 
is directly related to the occupation number of electrons 
calculated from the volume of the occupied states in the 
k space. 

Appendix B: Cauchy integral of using spline in- 
terpolation 

The Cauchy integral appears in various places in the 
DMFT calculation. Therefore, an efficient and accurate 
numerical calculation of the Cauchy integral is needed. 
We briefly explain a method of using the spline interpo- 
lation for the DOS. Let us assume that numerical data 
of DOS {yi} at energy points {x t } (i = 1, iVdatanumber) 
are given. In the cubic spline interpolation, 67 ' the DOS 
in the interval [x i7 x i+ i] is expressed as 

Ps,i( x ) = ^-{(^i+i - x) 3 M l + (x - Xi fM l+1 } 

Hy, - M^)^±L^ + {m+l _ M^±i)^i, (B .i) 



hi 



hi 



where hi = x i+ i — Xi. The quantity Mi is the second 
derivative of the DOS at a; = Xi, and is given by the 
usual procedure of the spline interpolation. 

The integral of the interval [x i} x i+ i] is calculated as 



rx i+1 

J Xi 



P ^x = (x l+1 -zf^- { z-x^ 
z — x 6 6 

+ {x t+1 -z)— (z-Xi) 



12 



12 



-p 3jJ (z)ln- 



Xi+\ 



Z Xi 



(B-2) 



The total integral is given by summing the contribution 
from each interval. Equation (B-2) is expressed as a com- 



bination of powers of quantities (x,+i — z) and {z — X,); 
therefore it is applicable for z when \z — Xi+ £ i + 1 1 j s not 
much lager than hi. When z approaches the edge of the 
integration, i.e., z —> Xi or ajj+i, the singularity of the 
logarithm term is removed by the counter contribution 
of the neighboring [xi_i,Xi] or [x i+ i,x i+2 ] region. The 
round-off error due to the subtraction of logarithm terms 
of neighboring regions is not serious even when z is ex- 
tremely near the edge point, because the divergence of 
the logarithm is very weak. The integral (B-2) is expected 
to be 0{z - £i±£i±i)-i when \z - Si±£i±l\ > There- 
fore, mutual cancellation occurs among terms in (B-2) in 
this limit. We find that the integral is re-expressed by a 
compact form in this case 



j-Xi + l 
J Xi 



P3,i( x ) 



dx 



z — X 



h\Mi 



oo , 

J K x l+1 -z> V+l 3(i/ + 3)(i/ + l) 
hi , v , Vt+i hfMi+i 



) 



(B-3) 



K z-Xi V + l 3(u + 3)(u + l)' 

This equation gives a highly accurate estimation of the 
integral even when terms are truncated up to v <~ 10. 

In most cases, it is convenient to use the linear inter- 
polation scheme with a very fine mesh for the DOS, 

p h i(x)=mi{x ) + , (B-4) 

where rrii = Vi +^~ Vi . The Cauchy integral is expressed 
by the following forms 

Pi,%{x) z-x i+1 
dx = —rriihi — pi,i{z) In , (B-5) 



z Xi 



and for \z — Xi+ f' +1 3> hi, 



Xi+l 



Pl,i( X ) j 1 u hi WfVi+i mjhj 



z — X 



dx = -i V [(_5i_)"(2«±l - ™ 



+1 



(^TO + ^)]- (B-6) 
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